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Abstract 

We examine the long-time behaviour of A + B(static) — > reaction-diffusion sys- 
tems with initially separated species A and B. All of our analysis is carried out for 
arbitrary (positive) values of the diffusion constant Z?a of particles A and initial 
concentrations do and bo of A's and B's. We derive general formulae for the location 
of the reaction zone centre, the total reaction rate, and the concentration profile 
of species A outside the reaction zone. The general properties of the reaction zone 
are studied with a help of the scaling ansatz. Using the mean-field approximation 
we find the functional forms of 'tails' of the reaction rate R and the dependence of 
the width of the reaction zone on the external parameters of the system. We also 
study the change in the kinetics of the system with D-q > in the limit D-q — * 0. 
Our results are supported by numerical solutions of the mean-field reaction-diffusion 
equation. 
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1 Introduction 



Investigation of the reaction fronts formed in the reaction-diffusion processes 
of the type A + B — > with initially separated species A and B has attracted 
a lot of recent interest. This is due to the fact that not only the kinetics of such 
systems exhibits many surprising features [1-5], but they are also amenable 
to experimental studies [3,6,7]. 
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A standard way to treat the initially separated problem analytically is to solve 
the partial differential equations [f] 

dpA _ n d 2 PA R 

(1) 

dpB _ n ^ 2 Pb d 

with the initial state given by 

p A (x,t = 0) = a H(-x) , 

(2) 

p B {x,t = 0) = 6 -^(^) , 



where p A (x,t) and Pb(^, are the mean local concentrations of A's and B's, 
R is the macroscopic reaction rate, H(x) denotes the Heavyside step function, 
and ao, bo, D A and Db are some constants related to the initial concentrations 
of species A and B and their diffusion coefficients, respectively. 

Equations (1) present the macroscopic approach to the problem. However, it 
has not been established yet how to relate the macroscopic reaction rate R 
to the microscopic picture of the initially separated reaction-diffusion systems 
below or at the critical dimension d c = 2. Dimensional [8-10] and renormal- 
isation group analyses [11-13] lead to the conclusion that above d c one can 
adopt the mean-field approximation R = kpxPB, with k being a constant. 
For 2D systems one expects logarithmic corrections to the mean-field picture 
[9,10,13]. One dimensional systems are usually studied by examining micro- 
scopic models in which, upon contact, the members of a pair A-B react with 
some probability p [10-18]. The analytical form of R(x,t) was derived for 
D B = [14], and for D B > and \x\ -> oo [12]. 

There are, however, several techniques which enables one to derive a lot of 
useful information from (1) even for d < d c , i.e. when the explicit form of 
R remains unknown. They are concentrated on the asymptotic, long-time 
behaviour of the reaction-diffusion systems, and include the renormalisation 
group analysis [11-13,17], the scaling ansatz [1,10], the quasistationary ap- 
proximation [8,19], and our approach of Ref. [20]. 

According to the scaling ansatz [1], the long-time behaviour of the reaction- 
diffusion system inside the reaction layer can be described with a help of some 
scaling functions Sa, Sb and 5r through 



p A (x,t) = VA r^S A (^-^) , (3) 
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PB (x,t) = VB r^s B (^-^ , (4) 

R(x,t) = VR t^S R (^^ ) j , (5) 

where Xf(t) oc t 1 ! 2 denotes the point at which the reaction rate R attains its 
maximal value, w(t) oc t a <C t 1 ^ 2 is the width of the reaction zone, rjA, ?7b 
and t^r are some parameters independent of x and t, and exponents a, f3, 7a 
and 7b are some positive constants given, for R oc paPb and nonzero diffusion 
constants, by a — |, (3 — | and 7a = 7b = §■ 

The quasistationary approximation [8,19] consists in the assumption that at 
sufficiently long times the kinetics of the front is governed by two characteristic 
time scales. One of them, tj oc (dlogj/dt) -1 oc t, controls the rate of change 
in the diffusive current J of particles arriving at the reaction layer. The other 
one, tf oc w 2 oc t 2a , is the equilibration time of the reaction front. If a < 1/2 
then t f /tj — > as t — > oo. Therefore, as t — > oo, in the vicinity of Xf the left 
sides of (1) become negligibly small compared to other terms. Consequently, 
if Da and D B are both nonzero, the asymptotic form of pa and p B inside the 
reaction layer is governed by much simpler equations 



(6) 



Da 1^ = R > 

which are to be solved with the boundary conditions 

dp^/dx — > — Ja(£), Pb —> as x — > — oo , 

(7) 

dp B /dx — > Jb(0) Pa - >• as rr — > +co . 



The most important feature of the quasistationary equations (6) is that they 
depend only on x, with time t being a parameter entering their solutions 
p\(x, t) and Pb(^, t) only through the time dependent boundary currents J a = 
Jb whose dependence on t, Da, D b , a and b has recently been derived 
analytically [20]. 

In a recent paper [20] we employed the quasistatic approximation to develop a 
new method of investigating the asymptotic kinetics of the initially separated 
reaction-diffusion systems. A peculiar feature of that approach is that it is 
concentrated mainly on the properties of the system outside the reaction zone. 
In this way, without imposing any special restrictions on the form of R, it 
relates exactly many quantities of physical interest to the values of external 
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parameters do, bo, D A and D B , which will enable us to investigate the limit 
D B —* analytically. In particular it was shown that there exist two limits, 
Cf = lim t _ >0O Xf(i)/-\/^ and Cj = lim^oo J{t)\fi. Given a , &o> -Da and D B , the 
value of Cf can by computed by solving 



$1 



C 



(8) 



where 



<3>(x) = [1 - erf(x)] exp(x 2 ) , 



(9) 



and erf(x) = 2ir 1//2 / x exp(— r/ 2 ) dr/ is the error function [21]. Then Cj can 
be calculated by solving 

C f = 2/Dlerf- 1 [(a - C A )/C A ] = 2 V^erf" 1 [(C B - b )/C B ] (10) 



and 



= C B W v exp -— - 



;n) 



where Ca and Cb are some constants controlling the form of pa and pb outside 
the reaction zone; specifically, for x <C Xf — w we have 



p A (x,t) = a -C A 



erf(x/yj4D A t) + 1 



;i2) 



and for x ^> Xf + w there is 



p B (x,t) = b + C B 



eri^x/y AD B t ] - 1 



(13) 



It was confirmed by several methods, including the renormalisation group 
analysis [11,12], numerical simulations [22] and heuristic arguments [20], that 
the values of the exponents a, (3, 7a and 7b as well as the form of the scaling 
functions S A , 5b an d 5r do not depend on D A , D B , a and b if the values 
of these parameters are nonzero. However, when one of the components is 
immobile (or 'static'), the asymptotic kinetics of the reaction front can change 
dramatically. For example, in the mean-field approximation the width of the 
front converges to a stationary value, and the reaction rate at Xf decreases as 
t -1 / 2 , which corresponds to a = and P — \ [4,22]. We have therefore two 
asymptotic universality classes: one characteristic for the 'dynamic' systems 
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in which both components diffuse, and the 'quasistatic' one observed if one of 
the diffusion constants is zero. Henceforth we shall assume Da > 0, so the two 
asymptotic universality classes will be distinguished by determining whether 
D B = or D B > 0. 

Although the peculiar kinetics of the systems with D B = was noticed quite 
early, so far nearly all of the research has been concentrated on the systems 
in which both of the diffusion constants Da and D B are nonzero. Closer in- 
spection of the powerful methods developed in the last decade to investigate 
the reaction-diffusion systems reveals that only the scaling ansatz, which is 
a purely mathematical concept, can be trivially employed to investigate both 
kinds of the systems. However, since even this relatively simple approach has 
so far been carried out only for the 'dynamic' problem [1,10], the theories of 
the two asymptotic universality classes remain practically disconnected. 

Several factors have lead to this situation. Theoretical results [14] derived for 
the case D B = utilized the basic property of such systems, immobility of 
particles B, in a way that cannot be extended for systems with D B > 0. On 
the other hand, most of the fundamental techniques developed for the case 
D B > has been based on the quasistationary approximation which requires 
that the ratio Ja/Jb of two opposite currents of particles A and B entering 
the reaction zone should asymptotically go to 1. This condition cannot be met 
by the systems with D B = 0, as in this case J B (t) = 0. 

The aim of our paper is to unify our understanding of these two kinds of the 
reaction-diffusion systems. The procedure we are using consists in detailed 
examination of the case D B = and comparison of the results with those 
already derived for D B > 0. We study the systems with D B = by means 
of various methods and for arbitrary (positive) values of a , b , and Da- In 
particular we show the counterpart of the quasistationary approximation (6) 
which should be used if D B = 0. We argue that the form of these equations, 
as well as their boundary conditions, determine the properties of the asymp- 
totic universality classes. In subsequent sections we study these properties for 
the systems with D B = using the heuristic theory of Ref. [20], the scaling 
ansatz, the mean-field approximation, and numerical analysis. Although some 
aspects of the systems with D B = turn out to be the same as of those with 
D B > 0, our analysis implies that these two cases should always be considered 
separately. 



2 The limit D B -> 

Consider equations (8) - (13) with the values of Da, cl and b fixed at some 
positive values, and D B going to 0. For physical reasons we expect that as D B 
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goes to 0, Cf converges to some nonzero limiting value, and so the argument 
of <!> on the right hand side of (8) diverges to infinity. We can therefore use an 
asymptotic property of the error function, x(l — erf(x)) exp(x 2 ) — > l/y/n as 
x — > oo [21], and reduce (8) to 



-C f 



v^oCf 



(14) 



Since $(x) diminishes monotonically from oo to as x grows from — oo to 
oo, the above equation has a unique, positive solution. Consequently, (10) and 
(11) imply that Cj and Ca also converge to some positive values, but Cb 
rapidly diverges to infinity. 

One can now use (10), (11) and (14) to arrive, after some algebra, at an 
important relation valid only if D B = 

6 C f = 2Cj. (15) 



This equation has a very natural physical interpretation. On the one hand the 
total number M(t) of reactions occurring by time t is asymptotically equal 

t t 
Jj a (t)&t~ JCj/y/r&T ~ 2CjVt. (16) 





On the other hand, however, for x ^> Xf + w we have p B ~ &o, and for 
iCif-towe expect pe ~ 0. Neglecting terms of order b w -C bot 1 ^ 2 we thus 
conclude that M(t) can as well be estimated by 

oo 

J [pbOM) - Pb(x, 0)] dx ~ b Xf(t) ~ b C{Vi, (17) 



which leads to (15). 

Our theory is consistent with the numerical simulations of Larralde et al [14], 
who considered the one dimensional system (i.e. for d < d c ) with D B = 0, 
D A = 1/2 and a /b = 1. They found that, respectively, for t = 500, 1000 
and 5000 the value of X{(t) was approximately equal to 11 ± 0.5, 16 ± 0.5 
and 36 ± 0.5, so that x f /\/i w 0.492 ± 0.022, 0.506 ± 0.016 and 0.509 ± 0.07, 
respectively. This is in excellent agreement with our equation (14), whose 
numerical solution reads Cf ~ 0.5060. Below, in Section 5, we will also verify 
our theory using the mean-field approximation (i.e. for d > d c ). 
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3 General consequences of the scaling ansatz 



Let D B = and a and b take on arbitrary (positive) values. Assume that 
the asymptotic solutions of the Galfi and Racz problem (1) in the long-time 
limit take on the scaling form (3) - (5) with a < 1/2. Inserting (3) and (4) into 
(1) and taking the limit t — > oo we find that at any x such that \x — Xf\ <C t 1 ! 2 
there is 

— ^— OC and D J^± oc r 7A-2a_ (lg) 

dt ox 



Hence, because a < 1/2, in the limit t — > oo the term dp\/dt becomes neg- 
ligibly small compared to D A d 2 p A j dx 2 . This implies that for \x — x^\ <C t 1 ! 2 
the form of the scaling functions can be determined by solving 

^ (19) 



The appropriate boundary conditions for these equations read 
dpx/dx — > —Jx(t), Pb — ► as x — > — oo , 

(20) 

Pa — >■ 0, p B — > &o as £ — ► +°o , 



where Ja(0 oc t l ^ 2 [20]. Equations (19) are the general counterparts of the 
quasistatic approximation (6) if Da > and Db = 0. 

The boundary conditions (20) determine the form of the boundary conditions 
for S A and S B except for a constant multiplier. We will take advantage of 
the fact that we are at liberty to multiply Sa and Sb by arbitrary constants 
(which can be compensated for by appropriate changes in f] A and t/b) and 
assume that the boundary conditions for Sa and Sb read 

S' A (z) — > —1, Sb(z) ^ as z — > — oo , 

(21) 

Sa(z) — > 0, Sb(z) — > 1 as z — > +oo , 



where the prime denotes the derivative with respect to z. Equations (21) 
immediately imply that 
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7b = 0, 



(22) 



VB = b . (23) 

The diffusive current J a of particles A for —y/Dpj C x C Xf - w is asymp- 
totically expected to be equal to Cj/y/i [20]. On the other hand, however, we 
can calculate it by inserting (3) into J a = —D\dp\/dx, which leads to 



7A + « = l/2, (24) 
r] A = C w Cj/D A , (25) 

where we denoted C w = \\m t ^ 00 w(t)/t a . 

Upon inserting the scaling ansatz into the first of equations (19) we come to 
7a + 2a = (3. Combining it with (24) we arrive at 

(3 -a = 1/2, (26) 



We thus see that the scaling ansatz imposes on the values of the scaling ex- 
ponents three relations (22), (24) and (26). Only the first of them takes on a 
different form if Db > 0, as in that case, by symmetry, 7a = 7b [10]. 

Equations (19) imply that inside the reaction zone 



Inserting into it the scaling ansatz (3) and (4), and carrying out our 'standard' 
limiting procedure (t — > oo at any x such that \x — Xf(t)\ <C t 1 ^ 2 ) we conclude, 
after some algebra involving (15), that Sa and Sb are related by a simple for- 
mula S'a(z) = S B (z). Upon integrating this equation and using the boundary 
conditions (21) to determine the integration constant we finally come to 

S> A (z) = Sb(z) - 1 . (28) 



Note that for D-q > the scaling functions Sa and Sb satisfy an entirely dif- 
ferent relation S\(z) = Sb(—z), which comes from the asymptotic symmetry 
of the system. Note also that relations (27) and (28) are independent of R 
and, consequently, of d. 



As a quick application of (28) consider a one dimensional system with Db = 0. 
It is known [14] that in this case S B = §(1 + erf(z)). So (28) and (21) lead to 



Sa(z) = |ierfc(z) = \ zerf(z) — z + tt exp(— z 2 ) 



(29) 
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4 The mean-field systems with D B = 



Consider a system governed by (1) with Da > 0, Db = and R = kp^PB- 
This particular form of R immediately implies 7a + 7b = P- Using now (22), 
(24) and (26) we conclude that 

a = 0, (3= 1/2, 7a = 1/2, and 7 B = . (30) 



These values are consistent with numerical simulations [4] and heuristic argu- 
ments [22]. 

It is easy to see that if length and time are measured in units of A = (-Da/ kbo) 1 ^ 2 
and r = 1/kao, respectively, then the solutions of (19) reduce to those ob- 
tained for the particular case D A — a — b — k — 1. Thus, in investigating 
the mean-field reaction-diffusion systems, it is sufficient to examine in detail 
only a system with some convenient values of the material parameters Da, 
ao, bo and k. The solutions for arbitrary values of these parameters can be 
then easily found by an appropriate choice of the units. This property guar- 
antees that any asymptotic length I satisfies /(-Da, a , b , k) = A/(l, 1,1,1). In 
particular, the asymptotic width of the reaction zone is given by 

W = lV ^ (31) 



where w denotes the asymptotic width of the reaction zone in the system 
with Da = a = b = k = 1. Numerical estimation of this parameter yields 
w ks 1.47 (see the next section for more details). 

Essentially the same line of reasoning was used in Ref. [20] to show that in 
the mean-field system with Db > the asymptotic width of the reaction 
zone is given by w = w 1 

(D> A D> B ) 1/3 (£:Cj)- 1 / 3 t 1 / 6 . However, in that paper it 
was incorrectly assumed that W\ = 1. We estimated the correct value of w\ 
numerically, obtaining w 1 1.38. 

Inserting now the scaling ansatz into (19) and using R = kpAPB we conclude 
that Sa and Sb satisfy 



S' a = w 2 SaSb , 
S b = w 2 S a Sb, 



(32) 
(33) 
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Upon inserting (28) into (32) we arrive at the nonlinear differential equation 
for the mean-field scaling function S\ 

S A = w 2 S A (S' A + 1) . (34) 

We can use it to estimate the behaviour of Sa(z) and Sr(z) for z ^> 1. In this 
region we expect S' A <C 1 and 5b ~ 1, so (34) reduces to S A = w 2 Sa, which 
implies 

Sn(z) oc Sa(z) oc exp(— wz) . (35) 

We can also investigate the tail of particles B which forms for z <C — 1. In this 
region we can assume Sa(z) ~ —z, so that (33) reduces to S' B = —w 2 zSb, 
which leads to 

Sb(z) oc exp(— \{wz) 2 ) and Sr(z) oc \z\ exp(— \{wz) 2 ) . (36) 

Thus, if Db = 0, the mean-field form of the scaling function S-r(z) is asym- 
metric, whereas for Db > this function is always symmetric, which is most 
easily seen in the symmetric case D A = D B and a = b . 

We will now investigate some properties of the limit D B — > which could not 
be analysed within the framework of the general theory presented in Section 2. 
As we already showed in our previous paper [20], if Db > then the mean- 
field density of particles B at X{ is asymptotically proportional to D B 2 ^ 3 t -1 / 3 . 
As this quantity has to be less than 6 0) we conclude that D B " 2//3 t~ 1 / 3 < const i. 
Therefore the time t* at which the mean-field system enters the long-time 
regime must satisfy 

t* > (£> B r 2 • const 2 . (37) 

Only for times satisfying this relation can one use the quasistatic approxi- 
mation (6). However, as D B — > 0, the right hand side of (37) diverges to 
infinity. Consequently, as Db goes to 0, t* diverges to infinity and in the lim- 
iting case Db = the kinetics of the system can never be described with the 
quasistatic approximation equations (6). Although we derived these conclu- 
sions only within the mean-field approximation, it is reasonable to expect that 
t* — > oo as Db — > in any initially separated reaction-diffusion system. 

To summarize the differences between the two asymptotic universality classes, 
in Table 1 we list their main properties in the mean-field approximation. The 
data for the case D B > come from References [1], [20] and [23]. 
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Table 1 



Comparison of the two asymptotic universality classes in the mean-field approxi- 
mation, 





D A > 0, Z?b > 


D A > 0, Db = 


Diff. eqs. for S A and Sb 

Diff. equation for 5a 

Sr(z), z > w 

St> I 7 1 7 <<^ ? / 1 

w(t) 


c/' _ c. c< 
D A — OA<->B 

S'l = S A (S A + z) 

z 3/4 exp( _| z 3/ 2) 

w^DaDb/Wj^H 1 / 6 


i^a — ^ uA'JB 

5 B = w 2 S a Sb 
S'l = w 2 S A (S' A + 1) 
exp(— wz) 
z exnf — i (wz)^") 

Wy/D A /kb 


a 


1/6 





(3 


2/3 


1/2 


7a 


1/3 


1/2 


7B 


1/3 






5 Numerical results 



To check the theory presented in the previous sections for the case D B = we 
solved numerically, using the finite-difference FTCS (Forward Time Centred 
Space) method, partial differential equations (1) with the mean-field reaction 
rate R = kp A pB- We present the data obtained for a® = 0.1, bo = 0.1, D A = 
0.1, and k = 0.02. Other values of these parameters yielded similar results. 

First of all we verified the theory presented in section 2. In Fig. 1 we show 
the plot of pa and pb in the vicinity of Xf at t — 10 7 . The dotted line was 
computed from (12). It perfectly matches the numerical solutions up to x ~ 
700, i.e. outside the reaction layer. Actually, in the region —1000 < x < 697, 
the relative error is less than 7-10~ 4 . Also, the value of C f computed from (14) 
is 0.2263, whereas its numerical estimation Xf/y/i for t = 10 7 is 0.2253. 

Next we investigated the scaling properties of the considered system. As some 
aspects of this problem were already considered [4], we present briefly only 
those results which are relevant to our paper. In Fig. 2 we show the log-log 
plot of w(t). We can see that it initially grows as y/i, which is a typical short 
time limit behaviour [24], but beginning from t ~ 10 3 it quickly converges to 
a constant value. This enabled us to estimate w ~ 1.47. 

Finally, in Fig. 3 we present the scaling plot of \/iR oc S R as a function of 
z = (x — Xf)/w, where we used (31) with w = 1.47. The plots for t = 10 6 
and 10 7 are practically indistinguishable. Note also two facts. First, for z > 
2 the semilog plot of y/tR is linear in z, in accordance with (35). Second, 
<Sr(z) is discontinuous at z — —Xf/w, which reflects the fact that Pb(x, t) is 
discontinuous at x — 0. 
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680 690 700 710 720 730 

X 



Fig. 1. The concentrations pA and pb of A's and B's in the vicinity of Xf pa 713 at 
t = 10 7 . The dotted line was computed from (12). 




1E+1 1E+2 1E+3 1E+4 1E+5 1E+6 1E+7 



t 



Fig. 2. The log-log plot of the width of the reaction front w as a function of time. 
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6 Summary and conclusions 



We have investigated the long-time behaviour of the concentrations pa and 
Pb of species A and B in the initially separated diffusion limited systems with 
Db = 0. All of our analysis was carried out for arbitrary values of Da, do and 
6o- 

First we derived the formulae (14) and (15) which together with (12) describe 
the behaviour of p\(x,t) outside the reaction zone. An interesting feature of 
these equations is that we expect them to be valid for any reaction-diffusion 
system of the type A + B (static) — > that exhibits the scaling behaviour with 
a < \. Thus we conclude that many important properties of the reaction- 
diffusion systems with D B = depend only on the values of a , b and D\, 
but not on the explicit form of R. This includes the location of the reaction 
zone centre (controlled by Cf), the total reaction rate (controlled by Cj), and 
the concentration profile of particles A outside the reaction zone (controlled 
by C\ and Cf). A similar situation is observed in the systems with D B > 
[20]. 

Next we investigated general consequences of the scaling ansatz. We concluded 
that it determines the value of 7b = and imposes relations (24) and (26) on 
the values of a, f3 and 7a- These relations are also valid for D B > 0. We proved 
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that the scaling functions Sa and Sb are related by a simple formula (28). We 
also concluded that for = the asymptotic forms of pa and pb inside the 
reaction front can be derived from equations (19) which are to be solved with 
time depending boundary conditions (20). Thus these equations determine the 
asymptotic properties of the reaction front. 



We also examined in detail the properties of the reaction zone in the mean- 
field approximation. In particular we determined the functional forms of Sa, 
Sb and Sr far from Xf and the dependence of the width of the reaction zone 
on k, do, bo and Da- 



Out analysis showed that the main differences in the behaviour of the initially 
separated reaction-diffusion systems with Db = and Db > arise from the 
fact that the term dpB/dt can be neglected in (1) only if the corresponding 
diffusion constant .Db is nonzero. Therefore, depending on whether Db is zero 
or not, the long-time behaviour of the reaction front is governed by entirely 
different partial differential equations. If Db > 0, then we use the usual qua- 
sistationary equations (6), otherwise we must employ (19). The different forms 
of these equations and their boundary conditions imply the different forms of 
their solutions and, consequently, different asymptotic properties of the two 
universality classes. Therefore the cases D B = and D B > should always be 
considered separately. 



There is, however, evidence that in some special cases the two asymptotic 
universality classes may be very much alike. Such surprising conclusion follows 
from the extensive numerical simulations of the one dimensional system carried 
recently out by Cornell [15]. He considered a system with D A = D B > and 
concluded that a = 1/4, f3 = 3/4, and that asymptotically R is a Gaussian 
centred at X{ with its width growing as t a . The same results were derived 
analytically for the one dimensional system with Db = by Larralde et al [14]. 
Note also that this form of R, or more generally — any R that depends on x 
and t explicitly rather than through pa{x, t) and Pb(x, t), uniquely determines 
the form of pA(x,t). This comes from the fact that whether or not Db = 0, 
asymptotical forms of pa and R are related by the same differential equation 
DAd 2 pA(x,t)/dx 2 = R(x,t) with the same boundary conditions. Therefore it 
is quite possible that the only difference between the one dimensional systems 
with Db = and D B > lies in the form of S B and the value of 7 B . It should 
be stressed, however, that the asymptotic kinetics of one dimensional systems 
with Db > has recently become the subject of controversial discussions [15- 
17,25,26], and further exploration of this topic is still required before the final 
conclusions can be made. 



14 



Acknowledgments 

This work was supported by University of Wroclaw Grant No 2115/W/IFT/95. 
References 

[I] L. Gain and Z. Racz, Phys. Rev. A 38 (1988) 3151. 

[2] Proceedings of the NIH Meeting on Models of Non-Classical Reaction Rates, 
J. Stat. Phys. 65, No. 5/6, 1991. 

[3] H. Taitelbaum, Y. E. L. Koo, S. Havlin, R. Kopelman and G. H. Weiss, Phys. 
Rev. A 46 (1992) 2151. 

[4] S. Havlin, M. Araujo, Y. Lereach, H. Larralde, A. Shehter, H. E. Stanley, 
P. Trunfio and B. Vilensky, Physica A 221 (1995) 1. 

[5] Z. Koza and H. Taitelbaum, Phys. Rev. E 54 (1996) R1040. 

[6] Y. E. Koo and R. Kopelman, J. Stat. Phys. 65 (1991) 893. 

[7] H. Taitelbaum in Diffusion Processes: Experiment, Theory and Simulations, 
Andrzej P§kalski (Ed.), Springer- Verlag, Berlin, 1994. 

[8] S. Cornell and M. Droz, Phys. Rev. Lett. 70 (1993) 3824. 

[9] P. L. Krapivsky, Phys. Rev. E 51 (1995) 4774. 

[10] S. Cornell, M. Droz and B. Chopard, Phys. Rev. A 44 (1991) 4826. 

[II] B. P. Lee and J. Cardy, Phys. Rev. E 50 (1994) R3287. 

[12] M. Howard and J. Cardy, J. Phys. A: Math. Gen. 28 (1995) 3599. 
[13] B. P. Lee and J. Cardy, J. Stat. Phys. 80 (1995) 971. 

[14] H. Larralde, M. Araujo, S. Havlin and H. E. Stanley, Phys. Rev. A 46 (1992) 
R6121. 

[15] S. J. Cornell, Phys. Rev. E 51 (1995) 4055. 

[16] M. Araujo, H. Larralde, S. Havlin and H. E. Stanley, Phys. Rev. Lett. 71 (1993) 
3592. 

[17] G. T. Barkema, M. J. Howard and J. L. Cardy Phys. Rev. E 53 (1996) R2017. 

[18] M. Hoyuelos, H. O. Martin and E. V. Albano, J. Phys. A: Math. Gen. 28 (1995) 
L483. 

[19] E. Ben-Nairn and S. Redner, J. Phys. A: Math. Gen. 28 (1992) L575. 
[20] Z. Koza, J. Stat. Phys. 85 (1996) 179. 

[21] Y. L. Luke, Integrals of Bessel Functions, McGraw-Hill, New York, 1962. 



15 



[22] Z. Jiang and C. Ebner, Phys. Rev. A 42 (1990) 7483. 

[23] H. Larralde, M. Araujo, S. Havlin and H. E. Stanley, Phys. Rev. A 46 (1992) 
855. 

[24] H. Taitelbaum, S. Havlin, J. E. Kiefer, B. Trus and G. H. Weiss, J. Stat. Phys. 
65 (1991) 873. 

[25] S. J. Cornell, Phys. Rev. Lett. 75 (1995) 2250. 

[26] M. Araujo, H. Larralde, S. Havlin and H. E. Stanley, Phys. Rev. Lett. 75 (1995) 
2251. 



16 



